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ABSTRACT 

We exploit the gravitational potential of a massive, rich cluster at z = 0.9 to study 
the internal properties of a gravitationally lensed galaxy at z=4.88. Using high res- 
■ olution HST imaging together with optical (VIMOS) and near-infrared (SINFONI) 

\ Integral Field Spectroscopy we have studied the rest-frame UV and optical properties 

of the lensed galaxy seen through the cluster RCS0224-002. Using a detailed gravita- 
tional lens model of the cluster we reconstruct the source-frame morphology on 200pc 
scales and find an ~L* Lyman-break galaxy with an intrinsic size of only 2.0x0.8 kpc, 



Or 



a velocity gradient of g60kms 1 and an implied dynamical mass of l.Ox 10 10 M G 



© 

within 2 kpc. We infer an integrated star-formation rate of just 12±2M yr" 1 from 
the intrinsic [On]A3727 emission line flux. The Lya emission appears redshifted by 
+200±40kms _1 with respect to the [Oil] emission. The Lya is also significantly more 
extended than the nebular emission, extending over 11.9x2.4kpc. Over this area, the 
Lya centroid varies by less than lOkrns -1 . We model the asymmetric Lya emission 
with an underlying Gaussian profile with an absorber in the blue wing and find that 
rN \ the underlying Lya emission line centroid is in excellent agreement with the [On] emis- 

■ sion line redshift. By examining the spatially resolved structure of the [Oil] and Lya 

emission lines we investigate the nature of this system. The model for local starburst 
galaxies suggested by Mass-Hesse et al. (2003) provides a good description of our data, 
and suggests that the galaxy is surrounded by a galactic-scale bi-polar outflow which 
has recently burst out of the system. The outflow, which appears to be currently lo- 
cated >30kpc from the galaxy, is escaping at a speed of upto ^500 kms -1 . Although 
the mass of the outflow is uncertain, the geometry and velocity of the outflow suggests 
that the ejected material is travelling far faster than escape velocity and will travel 
more than 1 Mpc (comoving) before eventually stalling. 

Key words: galaxies: high redshift; galaxies: starburst; gravitational lensing; galaxy 
clusters; Integral Field Spectroscopy; Gravitational Arcs: Individual: RCS0224-002 



1 INTRODUCTION 

One of the most important observational breakthroughs in 
recent years was the discovery that a significant fraction 
of hi gh-redshift galaxies are surrounded by "superwinds" 
(e.g. iPettini et al.1 1200I IShaplev etal] 120031 ; iBower et all 



120041 : 1 Wilman et al J 12003 ) - starburst and/or AGN driven 



outflows which expel gas from the galaxy potential, hence 
playing no further role in the star-formation history of the 
galaxy. This phenomenon is beginning to be understood by 



theorists as the missing link in galaxy formation models 
which are otherwise unable to match the shape and nor- 
malisation of the luminosity function l|Benson et al.l 120031 : 
iBaugh et ai]|200Sj h These feedback processes may also offer 
natural explanation as to why only 10 % of baryons cool to 
form stars (the Co smic Cooling Crisis; IWhite fc Reeslll978l ; 
iBalogh et aill200ll ). 

However, important questions remain unanswered. Ev- 
idence for these superwinds is usually based on observations 
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Figure 1. Le/i: True colour ffi>T VT-band image of the core of the lensing cluster RCS0224-002 at z=0.78. The contours mark the 
high-redshift critical curves (curves of infinite magnification) from the gravitational lens model described in §3. We also overlay the field 
of view of the SINFONI IFU (shown by the white box) which was used to map the [Oll]A3727 emission. The cluster galaxies which we 
are able to spectroscopically identify are labelled CG1-6. Rl is the radial counter-image of the z=4.88 arc and R2 is the 2=1.05 radial 
arc from Sand et al. (2005). Serendipitous background galaxies are labelled Al and A2 (VIMOS) and 1-4 (SINFONI) (see Appendix |A1. 
Right: VR(I+Lya) colour image of the cluster core generated from the VIMOS IFU datacube. The inner and outer curves show the 
2=4.88 caustic and critical curves respectively. The center of the cluster (0,0) is at a =02:24:34.255 6 =-00:02:32.39 (J2000) and we have 
rotated and aligned the HST and VIMOS data such that North is up and East is left in both panels. 



which compare the nebular emission line properties with the 
rest-frame UV emission and absorptio n lines (such as Lya, 
Ha and UV ISM absorption lines; e.g. lErb et al.lr2003l ). Ve- 
locity offsets of several hundred km/s have been measured, 
suggestive of large scale outflows comparable to starburst 
driven winds often observed in low-redshift Ultra-L uminous 
Infra - Red Galaxies (ULIRGs) in the local Universe l|Martinl 
2005). However, the current data lack spatial information, 
which is vital if we are to understand if material escapes 
into the Inter- Galactic Medium (IGM) or whether the out- 
flowing material eventually stalls, fragments and drains back 
onto the galaxy, potentially disrupting the disk and causing 
further bursts of star-formation. 

The key to resolving these issues is to identify similar 
features in the spatially resolved spectra of distant, young 
galaxies. However, at the redshifts where these feedback pro- 
cesses are at their peak activity (above z ~2, i.e. when galax- 
ies were most rapidly forming stars), the immense luminos- 
ity distances mean that the emission line fluxes of these dis- 
tant galaxies are dramatically reduced. This problem is com- 
pounded by the fact that the nebular emission lines (such 
as [On] , [Om] and Ha) are redshifted into the near-infrared, 
where the sky background is an order of magnitude brighter 
than in the optical. Coupled with the fact that galaxies at 
these redshifts typically have disk scale lengths of only a 
few kpc, even with adaptive optics assisted observations on 
eight to ten meter telescopes, only a handful of independent 
pixels can be recovered from a distant galaxy. 

By necessity therefore, above z ~ 2 most detailed 
studies of individual galaxies have concentrated on the 
most luminous (and therefore usually very active and 



massive) systems dBower et al.l 120041 : IWilman et alj 120051 ; 
ISwinbank etafl 120051 ). with very few "nor mal" galaxies 
being studied in detail beyond z = 2 (|Pettini et al] 
200al200d:lTeplitz et al]|200ol ; lFo"rster Schreiber et al]|200fl 
Genzel et al.ll20061 ). One way to overcome this problem is to 
use the deep potential of massive galaxy clusters to boost the 
flux and sizes of images of distant galaxies which serendipi- 
tously lie behind them. This natural magnification provides 
the opportunity to study intrinsically faint high redshift 
galaxies with a spatial resol ution that cannot b e attained via 



conventional observations (ISmail et al] 1 19961; Franx et al. 


19971; iTeDlitz et al.l2000l; Ellis et alj|2001 


;ICaniDusano et al. 


200 ll ISmith et al.N2002l; ISwinbank et al. 


120031; iKneib et al. 



12004 ; ISwinbank et al.ll2009 ). 



The natural magnification caused by the gravitational 
lens allows us to spatially resolve the morphologies and in- 
ternal dynamics of galaxies in a level of detail far greater 
than otherwise possible. For the galaxy discussed in this pa- 
per, the amplification factor of sixteen corresponds to an 
increase of over three magnitudes and makes it possible to 
study the galaxy on independent spatial scales of ~200pc 
at z =5. In contrast, without a gravitational lens, 1" corre- 
sponds to 6.5 kpc at 2=4.88. 

In this paper we present a VLT/IFU study of the z = 
4.88 gravitationally lensed galaxy behind the core of the rich 
lensing cluster RCS0224-002. This galaxy cluster was iden- 
tified as having a sig nificant concentration of red galaxies in 
the inner regions bv lGladders et al.l |2002). Follow-up spec- 
troscopy confirmed a redshift for the cluster of 0.773± .0021 . 
The ground based imaging from iGladders etafl (|2002D 
shows several tangential arc-like structures, many of which 
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Figure 2. Collapsed, one-dimensional spectrum of the arc from 
our SINFONI IFU observations. The [On] emission line has a 
redshift of 4.8757±0.0005 and a line width of o=180±30 kms" 1 . 
The inset shows the renormalised [On] emission lines for the two 
regions (A&B) shown in Figure 4. We also show a sky-spectrum 
(below) for comparison (scaled in flux for clarity). 

resemble images of lensed, background galaxies. One of 
the most striking of these is a multiply imaged arc ap- 
proximately 15" to the North-East of the Brightest Clus- 
ter Galaxy (BCG) (se e Fig. [I]). Follow-up spectroscopy by 
iGladders etafl (|2002r ) yielded a redshift of z=4.88 from 
identification of a strong Lya emission line at 7148A. 

We concentrate on the dynamics and star-forming prop- 
erties of the 2=4.88 arc observed through the rest-frame 
UV spectra with the VIMOS IFU and through the nebu- 
lar ([Oll]AA3726.1,3728.8A) emission line doublet observed 
with the SINFONI IFU in the near-infrared. The IFU data 
provide a two-dimensional map of the galaxy's properties 
in sky co-ordinates. In order to interpret this map we must 
correct for the magnification and distortion caused by the 
lensing potential using a mass model for the cluster lens. 
This model is constrained by requiring it to account for the 
positions and redshifts of the gravitationally-lensed arcs in 
the cluster. The lensing model then allows us to determine 
the source plane morphology and geometry of the continuum 
and line emission we use. 

The structure of this paper is as follows. In §2, we 
present the data on which this paper is based: a combina- 
tion of Hubble Space Telescope (HST), ground-based imag- 
ing, and integral field spectroscopy obtained with the SIN- 
FONI and VIMOS IFUs on the VLT. In §3 we present the 
source-plane properties of the arc and the spatially resolved 
spectroscopy and in §4 and §5 we present our discussion 
and conclusions respectively. Through-out this paper we use 
the Vega magnitude system and assume a cosmology with 
Ho = 72kms~ 1 , fi = 0.3 and A = 0.7. 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 HST imaging 

HST WFPC2 J 814 - and V 55 5~band observations of the 
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Figure 3. Image-plane observations of the ,2=4.88 arc: HST con- 
tinuum (greyscale), Lya emission (red contours) and [On] emis- 
sion (white contours). The image shows that the Lya emission 
is much more extended than either the [Oil] or continuum mor- 
phology. The solid bar in the top right hand corner represents 
the 0.8" seeing for both the VIMOS and SINFONI IFU obser- 
vations. Inset: Continuum subtracted, narrow-band image around 
the redshifted [On]A3727 emission of the arc from our SINFONI 
IFU observations. The highest surface brightness components in 
both the optical image and the [On] emission line map are well 
matched, though the [On] emission line map also shows that the 
highly sheared component has more structure than evident in the 
optical imaging. Right: 

lensing cluster RCS0224-002 were obtained from the HST 
public archivcQ- The Isn and Vsss-band observations were 
10.5 and 8.4 ks respectively and the data were reduced using 
the standard STSDAS package in iraf. The resulting image 
(Fig. [1]) covers the brightest central cluster galaxies and the 
target arc at a resolution of 0.0996" /pixel. From the HST 
imaging the z=4.88 arc is red (Veo6 — isi4=1.7±0.1) with an 
average surface brightness of /i/ sl4 =25.0 and an integrated 
magnitude of /=22.2. 

1 PID: 9135; Obtained from the Multimission Archive at the 
Space Telescope Science Institute (MAST). STScI is operated by 
the Association of Universities for Research in Astronomy, Inc., 
under NASA contract NAS5-26555. Support for MAST for non- 
HST data is provided by the NASA Office of Space Science via 
grant NAG5-7584 and by other grants and contracts. 
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Figure 4. Source-plane observations of the 2=4.88 galaxy from the HST, VIMOS and SINFONI observations. Top Left: [On] emission 
line intensity of the galaxy (dark regions represent regions of highest intensity). Top Right: [On] emission line velocity structure of the 
galaxy which shows a maximum velocity shift of 60±20kms~ 1 along the long axis of the galaxy. Bottom Left: Reconstructed true colour 
HST VI image of the 2=4.88 arc. Inset: Reconstructed HST image after a smoothing scale of 0.8" has been applied to the sky-plane 
image. Bottom Right: Reconstructed Lya emission line map of the 2=4.88 arc. The boxes show regions from which the spectra in Figure 6 
were extracted. 



2.2 VIMOS Integral Field Spectroscopy 

The z=4.88 arc w as observed with the VIMOS IFU 
ijLeFevre et al.ll2003l ) for a total of 43.2ks (split into 16x2700 
second exposures) between 2004 December 16 and 2005 De- 
cember 12 in <;0.6" seeing and photometric conditions. We 
used the MR-orange grism which results in a field of view 
of 27"x27" at 0.67" per pixel and a spectra resolution of 
A/AA ~1100 at 7000A. To reduce the data we used the 
VIMOS ESOREX pipeline which extracts, wavelength cali- 
brates and flat-fields the data. The final datacubes cover a 
wavelength range of 5000-1 1000A and has a resolution of 
6.6A FWHM (measured from the width of the sky-lines at 
a wavelength of 7150A). In all future sections, line widths 
are deconvolved for instrumental resolution. Flux calibra- 
tion was achieved by using observations of ESO standard 
stars. These observations were taken either immediately be- 
fore, or immediately after the science observations and were 
reduced in an identical manner. Since there are no point 
sources in our datacubes, we measure the seeing from the 
standard stars in the collapsed datacubes and derive typical 
seeing measurement of ~0.8" at a typical airmass 1.1. Dur- 
ing the standard star observation, the ESO seeing monitor 
registered median seeing of 0.7" - hence we conservatively 
assume that the VIMOS IFU observations were taken in 0.8" 
seeing. 

After reducing each of the observations to a datacube, 
we construct and apply a master flat-field in order to cor- 
rect for the fringing effects which are particular prominent 
above 8000A. In each of the (spatially dithered) datacubes 



we construct a broad-band image by collapsing the datacube 
in wavelength between 6000A and 7500A and mask out the 
the brightest objects. We then apply this mask at each wave- 
length and stack the resulting datacube (without spatial off- 
sets and ignoring masked pixels). This produces a "master- 
flatfield" which is then divided into each of the science ob- 
servations to reduce the effects of the fringing. These indi- 
vidual datacubes are then sky-subtracted by masking the 
brightest objects in the individual frames and constructing 
a sky-spectrum by collapsing the datacube in both spatial 
dimensions (since the four quadrants of the VIMOS IFU 
have slightly different throughput's, we note that the sky- 
subtraction was performed quadrant-by-quadrant) . To build 
the final datacube, we create broad band images from each 
datacube and cross correlate these in order to spatially align 
and create the final mosaic. The final datacube is created us- 
ing an average with a 3a clip to reject cosmic rays. 

In Fig. [T] we show a colour image of the cluster and 
arc, generated from the datacube. We generate three colours 
from the datacube by collapsing the datacube between 5000- 
6000A, 6000-7000A and 7000-8000A as blue, green and red 
respectively. Alternatively, the data-cube can be viewed as 
a sequence of wavelength slices to make an animation (this 
can be viewed at: 

http://star-www.dur.ac.Uk/~ams/lensing/RCS0224_ z5/l ) 
2.3 SINFONI IFU Observations 

To map the nebular emission line properties of the 2=4.8 8 
arc, we used the SINFONI IFU jEisenhauer et al.ll2003h . 
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Figure 5. Collapsed, one-dimensional spectrum of the arc from 
our VIMOS IFU observations. The dashed vertical lines show the 
expected position of the UV emission and absorption lines for 
the systemic redshift of z =4.8757 (as measured from the [On] 
emission line). The Lya emission appears redshifted with respect 
to the systemic velocity, and has an asymmetric line profile. We 
also detect the weak UV-ISM lines of Siiv. These lie in a region 
devoid of strong sky emission. The hashed regions show regions 
of strong sky emission and the horizontal dashed line marks a 
continuum level of zero. Top Inset: The Lya emission line with 
the model described in §3.2.2. The black (solid) line shows the 
observed spectra. The asymmetric line is best fit with an under- 
lying Gaussian emission line profile (blue dotted line) combined 
with a Voigt profile absorber (green dashed line) and a third 
Gaussian which accounts for the red wing of emission. Whilst 
the observed Lya emission appears redshifted from the systemic 
velocity, the best-fit underlying emission profile has a centroid in 
excellent agreement with the [On] emission. Lower Inset: Spectra 
around the rest-frame UV absorption lines SilvAA1393.76,1402.8. 
We overlay the best-fit template (blue) from Shapley et al. (2003) 
and use this to infer a velocity shift of -400±100 km s~ 1 from the 
systemic redshift. 



At 2=4.88 the [On] A A3726. 1,3728.8 emission line doublet is 
redshifted to 2.19^im and into a stable part of the A"-band 
relatively free from strong OH airglow emission. The SIN- 
FONI IFU uses an image sheer and mirrors to reformat a 
field of 10x10" at a spatial resolution of 0.25"/pixel. We 
used the HK grating which results in a spectral resolution 
of A/AA=1700 at 2.20/im (the sky-lines have a resolution 
13A FWHM at this wavelength) and covers a wavelength 
range of 1.451 to 2.463/xm. In all future sections, emission 
line widths are deconvolved for the instrumental resolution. 
To observe the target we used ABBA chop sequences, how- 
ever since the arc only filled one half of the datacube, we 
chopped 2" West to sky (see Fig. [JJ , thus keeping the ob- 
ject inside the IFU. We observed the target for a total of 
43.2 ks (split into 16x2700 second exposures) between 2005 
August 10 and 2005 September 28 in £0.6" seeing and pho- 
tometric conditions. Individual exposures were reduced us- 
ing the SINFONI ESOREX data reduction pipeline which ex- 
tracts, flatfields, wavelength calibrates and forms the dat- 
acube. The final datacube was generated by aligning the 
individual data-cubes and then combining the using an aver- 



age with a 3<r clip to reject cosmic rays. For flux calibration, 
standard stars were observed each night during either im- 
mediately before or after the science exposures. These were 
reduced in an identical manner to the science observations. 
In Figure [2] we show the one-dimensional spectrum of the 
arc (generated by collapsing the datacube over the object 
in both spatial domains) and in Figure [3] we also show a 
narrow-band image of the arc generated by collapsing the 
datacube in the wavelength direction over the [On] emission 
line. Figure [3] also shows the HST image of the arc with the 
contours from the [On] and Lya overlaid. 



3 ANALYSIS & RESULTS 
3.1 Mass Modelling 

Reconstruction of the intrinsic properties of the galaxy at 
.z=4.88 requires removal of the gravitational magnification 
from the observables. We therefore use the four spectroscop- 
ically confirmed images of the 2=4.88 galaxy to constrain a 
model of the mass distribution in the cluster. These four im- 
ages c omprise the three ima ges in the tangential arc identi- 
fied bv lGladders et al.l ((2002) and the fourth image identified 
via our integral field spectroscopy (Appendix |A"| . The can - 
didate radial arc at 2=1.05 (A3 in Fig. IT1 ISandet al.ll2005h 
and the object identified at 2=3.66 adjacent to the 2=4.88 
arc (Al in Fig. [1]) appear not to be multiply-imaged, and 
so are not included as constraints. The configuration of the 
four observed images and the lack of strong radial amplifica- 
tion of the fourth image implies that the fifth image lies very 
close to the centre of the cluster mass distribution. We there- 
fore expect the fifth image to be strongly de-magnified and 
thus likely undetectable in the current data. We identify the 
dense knot of emission in each of the four observed images 
as being images of the same underlying region of the galaxy, 
and use the position of t hese knots to const rain the model. 
Following Appendix A of lSmith et al. I (|2005h . the number of 
constraints is therefore nc=6. 

The central region of RCS0224 contains two bright 
elliptical galaxies: CGI and CG2 (Fig. [TJ - we there- 
fore construct a three component model comprising the 
cluster-scale mass distribution (dark matter and gas), 
CGI and CG2. We parametrise all three mass com- 
ponents as truncated pseudo-isothermal elliptical mass 
distributions, each described b y the following parame- 
ters: {x, y, e, 0, rc o re , i"cut, «di SP | jKassiola fc Kovnerl Il993l : 
iKneib et all 1 19961 ) - i.e. 21 parameters in total, comprising 
the center, ellipticity, position angle, core radius and cut-off 
radius of each mass component. To obtain a well-constrained 
model, (and given the number of constraints above: nc=6) 
we prefer no more than 6 of these 21 parameters to be free 
parameters in the lens model. The geometrical parameters 
(x, y, e, 9) describing CGI and CG2 are therefore matched to 
their 2-dimensional light distributions, and their respective 
r corc and r cut parameters are matched to the shape of the 
light profile of each galaxy. Finally, the mass-to-light ratio 
of each galaxy is fixed by matching their velocity dispersions 
to the measured velocity dispersion (Table |B)| . Turning to 
the parameters describing the cluster-scale mass, we fix both 
r corc and r cut at 75 kpc and 1 Mpc respectively. The former 
value is representative of massive clusters, and our decision 
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to fix this parameter is caused by the non-detection of the 
central fifth image of the z=4.88 arc. The lack of the fifth 
image reduces our ability to constrain the radial shape of 
the mass distribution within the tangential critical curve. 
We also run models with r CO rc=50kpc and 100 kpc to con- 
firm that the final results are not altered by the adopted 
value of r coro . For completeness, we include the result of 
this check in the error on the lens magnification calculated 
below. Fixing r cut = l Mpc is motivated by the small field of 
view of the HST data and within those data the small an- 
gular scale subtended by the strong lensing constraints. We 
are therefore unable to constrain this parameter, but our 
results are insensitive to its value. In summary, we therefore 
fit a model with five free parameters to the data. The free 
and fixed parameters are listed in Table 1. 

All of the len s modelling is performed us ing the 
lenstool software (|Kneidfl993l ; iKneib et al.lll996l) . incor- 
porating a Markov Chain Monte Carlo (MCMC) sampler 
(Julio et al. 2006, in preparation). The model parameters 
with the lowest % 2 found with this MCMC approach and 
68% confidence intervals around these parameters values af- 
ter marginalising over the other four parameters in each 
case, are listed in Table 1, and we show the critical curve 
of the best fit model in Fig. [T] The lens model predicts the 
fifth image to lie at the center of CGI and to be very faint 
- ife36, i.e. de-magnified by ~14 magnitudes relative to the 
tangential arc. Although the model makes a clear prediction 
for the location of this image, it is clearly too faint to de- 
tect with the available data. We ray-trace the z=4.88 galaxy 
through the family of models within the 68% confidence sur- 
face in the 5-dimensional parameter space to compute the 
mean, luminosity weighted magnification of /i=16±2 (which 
corresponds to Am=3.0±0.2 magnitudes). 

Accounting for the lensing magnification, the intrinsic 
magnitude of this galaxy is 1=25.2. In comparison, deep 
imaging surveys of Lym an-break galaxies a t z~5 have de- 
rived an L* of i ~25.3 (jOuchi et all [2004). This suggests 
that the z =4.88 galaxy is typical of the UV continuum se- 
lected LBG population at this redshift. In Fig. [4] we show 
the reconstructed HST image of the galaxy. In the source- 
plane the galaxy has a FWHM of only 2.0x0.8 kpc, and in- 
dicates a flattened (or bar-like) geometry. Such an elongated 
morphology is not unusual for galaxies at these redshifts: in 
a recent high-resolution imaging survey of LBGs at z>2.5 
with HST, Navindranath et al. (2006) concluded that upto 
50% of LBGs show evidence for bar-like morphologies (with 
mean scale lengths of 1.7-2.0 kpc). Thus the galaxy appears 
typical in terms of its source plane morphology, brightness 
and size of the population at z ~ 5. 



3.2 Spatially Resolved Spectroscopy 

3.2.1 Nebular emission 

By fitting a double Gaussian profile of fixed separation (cor- 
responding to the redshifted separation of the doublet) and 
intensity but variable width, redshift and intensities of the 
doublet to the collapsed spectrum, the [On] emission yields 
a systemic redshift of 4.8757±0.0010 and and intrinsic width 
of <r=100±20 km s _1 (in the rest-frame of the galaxy). We 
assume that the [On] maps the systemic redshift and all 
velocities are given with respect to this value. The inte- 



grated emission line flux of 5±lxl0~ Wm~ suggests a 
star-formation rate (corrected for lensing amplification, but 
uncorrected for reddening) of 12±2M q yr~ (assumin g the 



calibration of L([On]) to SFR given in lKennicuttj|l998f ). 

In Figure [2] we show the intensity distribution and col- 
lapsed spectrum of the [On] emission. Intensity and velocity 
maps were derived by fitting the [On] emission line dou- 
blet in each 0.25" pixel using the same technique described 
above. We used a \ 2 minimisation procedure, taking into 
account the greater noise of at the positions of the sky-lines. 
In cases where the fit failed to detect the line the region was 
increased to 2x2 pixels (0.5" xO. 5"). Using a continuum fit, 
we required a minimum signal-to-noise of 3 to detect the 
line, and when this criterion is met, we fit the [On] emis- 
sion line doublet with a Gaussian profile of fixed separation 
allowing the central wavelength and normalisation to vary. 
In Figure [4] we show the intensity distribution and velocity 
structure of the [On] emission line flux. 

It is clear that the brightest components in the rest- 
frame UV are also the brightest in the [On] emission. In 
the source plane we find that the eastern most compo- 
nent is marginally redshifted relative to the western com- 
ponent. Relative to the combined spectrum, the eastern re- 
gion has v = —20 ± 20kms~ 1 (with a velocity dispersion 
of <T=80±18kms~ 1 ), while the western component is at 
+20±15 kms -1 (and a velocity dispersion of 85±12kms _1 ). 
The quoted uncertainty here is dominated by possible vari- 
ations in the [On] doublet line ratios, which was left free in 
these fits. 

Since the [On] velocity dispersion (a) reflects the dy- 
namics of the gas in the galaxies' potential well we estimate 
the dynami cal mass o f the g alaxy. Following the same pre- 
scription of Erb et al l l|2006tt , we use the velocity dispersion 
and spatial extent of the [On] flux to infer a dynamical mass 
of ~ 1 x 10 10 M within a radius of 2 kpc. Clearly this 
mass has large uncertainties since the mass depends on the 
mass density profile, velocity anisotropy and relative contri- 
butions to a from random motions or rotation and possible 
differences between the total mass and that of the tracer 
particles used to measure it. Nevertheless, it is worth noting 
that the mass we derive is a factor of ~ 5 x smaller than the 
median dynamical mass of z~3 LBGs from lErb etafl (|2006T ) 
measured in exactly the same way. 

3.2.2 Lya 

In Fig [4] we show the reconstructed source-plane Lya mor- 
phology of the arc. It is immediately apparent that in the 
source-plane, the Lya is more extended than the [On], or 
the continuum light, extending over 11.9x2.4 kpc (FWHM). 
Figure [S] shows the collapsed spectrum of the arc overlaid 
with the expected position of the Lya emission line for 
the systemic redshift (we also mark the expected position 
of the UV ISM lines which lie in regions of the sky free 
from strong OH airglow emission). The integrated Lya line 
flux is 3.7±0.5xl0" 16 ergs _1 cm 
formation rate of 5±2M© yr _1 

the 1500A co ntinuum flux sug gests a star-formation rate of 
3±2M yr- 1 l |Kennicutdli99af l. 

From the collapsed spectrum of the galaxy we mea- 
sure a redshift from the centroid of the Lya emission line 
and measure a velocity offset from the systemic redshift of 



2 , corresponding to a star- 
(unlensed). In comparison 
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Table 1. 

Gravitational Lens Model Parameters 



AR.A ADcc e & r COTC r cut ^diap 

(") (") (deg) (kpc) (kpc) (km/s) 



DM halo -0.6±1.6 2.1±1.4 0.18±0.02 -49±5 [75] [1000] 935±30 
CGI [0.0] [0.0] [0.07] [169] [0.2] [42] [230] 

CG2 [2.4] [-1-2] [0.11] [202] [0.2] [38] [220] 



Table 1. Note: Numbers in square brackets are fixed in the fit. Position angles are clockwise from North. 



+200±40 kms -1 . Velocity offsets of this magnitude between 
Lya and nebular emission lines (such as [Oh] or Ha) are a 
common featur e of galaxies at z ~3 l|Shaplev et all 120031; 
lErb et all I2003T ). In both the collapsed spectrum (Fig. [5j) 
and the spatially resolved spectra (Fig. O the strong Lya 
from the galaxy shows an asymmetric profile in which the 
blue wing of the Lya is absorbed, presumably giving rise to 
the apparent redshift offset between the [On] and Lya. It is 
also clear that there is an extended wing of emission in the 
red side of the line. 

Asymmetry in the Lya emission line spectra of LBGs 
at z ~2-3 are common and often give rise to velocity offsets 
of several hundred km/s. These velocity offsets are usually 
attributed to galactic-scale outflows, presumably powered 
by supernovae. However, since the optical depth of the Lya 
forest is an order of magnitude larger at z = 5 than at z = 3, 
(where most spectroscopic studies have taken place) we first 
investigate whether the IGM could cause the absorption in 
the blue wing of the emission line. 

At z = 5 the mean optical depth of the Lya forest is 
1.7±0.2 l|Fan et al.| [2006) which corresponds to a transmit- 
ted flux of 20%. We therefore investigate whether the Lya 
profile could be explained as the result of average IGM ab- 
sorption, or whether the absorption needed to be due to ma- 
terial associated with the lensed galaxy. The mean optical 
depth is measured by averaging the absorption of individual 
clouds on scales of 6000 km s -1 (6Mpc comoving). We fit 
the Lya emission line with an underlying Gaussian emission 
line profile (centred on the systemic redshift of the galaxy 
as measured from the [On] emission), convolved with a step 
function which models the average IGM (all profiles are then 
convolved with the instrumental resolution) . We also include 
a broad Gaussian emission line profile (also centred at the 
systemic redshift) to account for the extended red wing of 
emission. Allowing the redshift of the IGM component to 
vary we find that an acceptable match emission line shape 
is best matched when the IGM step function is placed be- 
tween -80 and +10kms _1 from the galaxy (between -85 kpc 
and +10 kpc comoving for out adopted cosmology). This is 
inconsistent with our assumption of a large-scale average ab- 
sorption, implying that the absorption we see results from an 
individual cloud (or clouds) lying close to the target galaxy. 

Moreover, a simple step function fails to adequately re- 
produce the line shape of the blue wing. A better fit is 
obtained if the line shape is instead modified by a single 
absorbing cloud. We model the emission profile by Gaus- 
sian profile combined with a Voigt profile absorber in the 
blue wing. In this fit the wavelength of the underlying 
best fit Gaussian profile (blue dotted line in Fig [5j) was 
allowed to vary, but the fitted centroid at 7143.0±0.5A 
(z =4.8758±0.0001) is in excellent agreement with the neb- 



ular emission line redshift as measured from [On]. This 
underlying Gaussian profile also has a width of cr=4.8A 
(cr=200kms _1 in the rest-frame of the galaxy). The ab- 
sorber in the blue wing of the Lya emission (green dashed 
line) has a centroid of 7138.0±0.lA (z=4.8717±0.0001, or 
Av = — 500kms _1 relative to the intrinsic emission) and 
a width of ~100kms -1 . Since the flux in the Voigt pro- 
file reflects the column density of the neutral gas, we infer 
N m = 1.6±f:5 x 10 19 cm~ 2 and cr=110±50 kins' 1 . The large 
uncertainties in the value of Nm and a reflect that fact that 
the inferred value of Am lies between the flat and linear part 
of the curve of growth for Hydrogen. 

With this two component model we are able to fit the 
asymmetry in both the collapsed and spatially resolved spec- 
tra. However, even with two components there remains ex- 
cess flux in the red wing of the emission line. To compensate 
for this emission, we also fit a third Gaussian profile compo- 
nent in the red wing, again allowing the redshift, width and 
intensity to vary. The improvement between the three fit 
and two fit models is A\ 2 — 12 — 38 (~ 3 — 6a; depending 
on whether we fit the individual spectra or the combined 
spectrum). The best fit model therefore includes a third 
component with redshift centroid at 7153.44 (z=4.8860) and 
a width of cr=260kms~ 1 . The velocity offsets of the three 
components are summarised in Table 2. As we discuss in 
§13] this third component may naturally arise from Lya pho- 
tons back-scatter on the interior of an outflowing shell. Thus 
the best-fit model to the asymmetric Lya emission is one in 
which the Lya photons are generated in similar regions (and 
at the same redshift) as the [On], but the emission profile 
is modified by foreground neutral material. This foreground 
material absorbs/scatters photons from the blue wing of the 
line, causing an apparent redshift offset between the nebular 
emission and Lya. 

The description of the line profile given above is typical 
of the emission profile seen in other Lyman break galaxies. 
The key advantage of observing the lensed system is that 
we are able to spatially resolve the emission line, and hence 
investigate the spatial variations in the absorption and emis- 
sion. We divide the IFU data cube into the spatial regions 
shown in Fig. [4] For reference these have been labelled 0-15. 
The resulting spectra are shown in Fig. [6] The most striking 
feature of this figure is that (apart from variations in inten- 
sity) the line profile appears remarkably constant across the 
galaxy. For each spatial element, we measure the centroid 
of the emission and find that the intrinsic emission line cen- 
troid varies by less than lOkrns -1 (rest-frame) across the 
whole galaxy image (this value was derived by both fitting 
the 2 and 3 component models above). Indeed, as Fig. [5] 
shows, the asymmetry in the Lya profile is evident across 
the whole galaxy. Applying a free fit, each of the individ- 
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Figure 6. One-dimensional spectra around the redshifted Lya 
emission from the 2=4.88 galaxy from the fifteen regions shown 
in Figure 4 (black). Each of the spectra has the composite spec- 
tra (scaled) and overlaid for comparison (red dashed). The Lya 
emission in all of the spectra show an asymmetric profile, with 
the blue wing of the Lya emission truncated making it appear 
redshifted by +200±40kms _1 from the nebular ([On]) emission 
line. Furthermore, the centroid of the Lya emission varies by less 
than lOkrns - 1 (rest) across the whole galaxy image. 

ual spectra is consistent with the mean, with the maximum 
Ax 2 =l-3 (and an average A% 2 ~0.5). From the position of 
the Voigt absorber, we also note that the blue-edge of the 
emission line cuts off at a wavelength corresponding to vari- 
ations of less than 20kms" 1 at a redshift of 4.871 7±0.0001. 
This constancy in the line shape and centroid across the 
galaxy image is discussed further below. 

3.2.3 UV absorption 

In the integrated rest-frame UV continuum spectrum of the 
galaxy, we also detect the weak Silv(AA1393,1402.77) ab- 
sorption features. Whilst these are very weak in our spec- 
tra (the S/N in continuum in regions free from strong sky 
emission is only ~3), the separation between the doublet 
is 9.80±0.15A (rest), is in excellent agreement with the 
expected value of 9.77. To calculate the redshift of these 
absorption lines we cross correlate the spectrum around 
8200A with the rest-frame UV composite spectrum of Ly- 



Table 2. 

Properties of the emission and absorption lines 



Component 


z 


Av 

km s — 1 


km s~ 1 


V shear 

kms -1 


[On] 


4.8757[5] 





100±30 


-30< v <30 


Lya (Principle) 


4.8760[5] 


+20±50 


230±40 


£20 


Lya (red wing) 


4.8860[f5] 


+520±150 


260±80 


£50 


Lya (blue abs. ) 


4.8680[10] 


-500±100 


110±30 


£20 


Siiv absorption 


4.8700[10] 


-400±100 


~480 


£150 



Table 2. Note the value given in the [] z column is the error 
in the last decimal place. v ahf , ar denotes the maximum velocity 
gradient/shear in the given emission/absorption line. 

man break galaxies from Sha plev et al.l l)2003h and find that 
these absorption lines are blue-shifted by — 400±100kms -1 
(2=4. 870±0. 0007) from the systemic redshift as defined by 
the [On] emission line (Figure [SJ. It is worth noting that 
in the LBG composite spectrum from lShaplev et"al] (|2003l ) 
the interstellar absorption features (OlA1303+SillA1260) are 
detected with similar significance as the Siiv. However, at 
z = 4.88 these are redshifted to 7656A and into the middle 
of the Frauenhoffer A-band and therefore are therefore not 
expected to be detected. 

The observed velocity offsets between the nebular 
([Oil]), Lya and UV ISM lines are therefore consistent with 
the situation seen in most high-redshift Lyman Break Galax- 
ies where veloci ty offsets of seve r al hundred km/s ha ve been 
measured (e.g. lErb et al.|[2003l ; ISteidel et al.ll2004) . These 
velocity offsets are attributed to starburst driven superwinds 
and we base much of our interpretation on this model. 



4 DISCUSSION 

Before attempting to interpret the observations, we briefly 
review the results of the observations. The [On] emission 
maps the star-forming regions, which appear to have an elon- 
gated (bar-like) morphology with a spatial extent of only 
2.0x0.8kpc, which appears to contain all of the stars in the 
galaxy. We use the [On] emission to infer a systemic red- 
shift of 4.8757±0.0005. Using the [On] emission we infer an 
integ rated star- formation rate of 12±2M0yr _1 ( Kc nnicuttl 
ll998T ) and the spatial extent and line width of the [On] sug- 
gest a dynamical mass of ~1 x 10 10 M© within 2 kpc. Turning 
to the rest-frame UV spectrum: the Lya emission is much 
more extended than the [On], extending over 11.9x2.4 kpc 
(i.e. at least a factor of 4x the spatial extent of the [Oil]). 
Moreover, the peak of the Lycv emission is redshifted by 
+200±40kms~ 1 from the [On]. Across the whole spatial 
extent, the Lya emission line centroid varies by less than 
10 km s" 1 (rest-frame) and the line profile is remarkably con- 
stant over the same spatial extent. The Lya emission (both 
in the collapsed spectrum and spatially resolved) has an 
asymmetric profile which is best explained by a three com- 
ponent model. The underlying Gaussian emission line profile 
has a centroid in excellent agreement with the nebular emis- 
sion, an absorber which is blue-shifted by — 500±100 kms -1 
from the systemic redshift and an extended "red wing" of 
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emission is redshifted by +520±150 kms -1 . Finally, the UV 
ISM lines of Siiv are blue-shifted by -400±150kms _1 from 
the systemic redshift, although it has not been possible to 
derive any spatial information across the galaxy from these 
weak lines. 



4.1 A starburst with a bi-conical outflow? 

The observed velocity offsets between the [On] , Lya and UV 
ISM lines in RCS0224arc are consistent with situation seen 
in most high redshift Lyman Break Galaxies in which the 
presence of "superwinds" (driven by the collective effects of 
star-formation and supernovae) naturally explain the obser- 
vations. Our spatially resolved spectra therefore offer us the 
opportunity to study the spatial structure and energetics of 
the outflowing material in this galaxy. 

The velocity offset between the [On] and transmit- 
ted Lya emission profile are usually interpreted in terms 
of a starburst driven outflow. This model provides a 
good description of the integrated line emission profiles 
of high redshift galaxies and more detailed observations 
of lo ca l star bursts (eg., Pettini et al.l 20021; Shaplev et al.l 



20031; iTenorio-Tagle et all Il999l ; iHeckman et all I2060F 
Grimes et al.ll2006l) . Herewe investigate the implication for 
this model from the spatially resolved profile that we have 
ob served. We base our int erpretation on the specific models 
of iMas-Hesse et all (|2003h which assume that the spatially 
compact [On] nebular emission traces the star forming re- 
gions in the nucleus of a star-bursting galaxy. 

In these models, the velocity offsets between [On] and 
Lya is generated by a star forming region embedded in an 
outflowing bubble. The models have been extensively inves- 
tigated in local starburst galaxies where similar spatially 
resolved data is available. Indeed, there is strong similar- 
ity between t he data we presen t , and the local galaxy Haro 
2 studied by IMas-Hesse et al] 1 2003h . The spatial extend 
(500 pc) and velocity shift (w s heli ~ 200 kms -1 ) seen in Haro 
2 are smaller than the z — 4.88 galaxy but otherwise the sys- 
tems are remarkably similar. Both show Lya emission that 
extends beyond the nebular emission region and both show 
a blue absorption cut-off to the Lya emission line and a sim- 
ilar red-wing. Crucially, both data sets show no detectable 
variation in the shape of the Lya emission profile across 
the emission region. It is remarkable that this galaxy, seen 
1 Gyr after the Big-Bang is so similar to the local system at 
z = 0.0049. 

The basis of the interpretation is that the star-burst 
generates a large over pressurised region which surrounds 
the nebula. Within the ISM, the flow is Reighly- Taylor un- 
stable so that the initial shell breaks up leaving behind a 
denser gas clumps that may be responsible for the extended 
Lya emission region. The outflow tends to break out of the 
galaxy along the axis of least resistance. In the case of a 
disk (or flattened) galaxy this generates an approximately 
bi-conical outflow. When the outflow reaches the smoother, 
lower density gas distribution in the halo of the galaxy, a 
new shell develops as the outflow sweeps up material from 
the galaxy's halo, creating a screen between the observer 
and the source. As the shell travels away from the galaxy, 
it becomes den se enough to recombine and its Lya opacity 
becomes large. IMas-Hesse et all (2003) interpret Haro 2 as 



being an example of a system in these later stages, and the 
z = 5 arc fits into the paradigm well. 

The model interprets the spatially integrated proper- 
ties of the emission and absorption lines seen in Figure [6] as 
follows: the main peak of Lya emission comes from photons 
emitted from the star forming regions. To reach the observer 
these must pass through part of the foreground (blue shifted) 
shell. This absorption causes the peak emission wavelength 
to appear redshifted relative to the nebular emission lines. 
In addition, photons may be scattered or emitted from the 
receding shell. Photons that are either created on the inner 
surface of the shell (e.g. by UV irradiation from the star- 
burst), or multiply scattered within the receding shell so 
that they acquire the mean ve locity of the shell (see the dis- 
cussion ofH anscn fc Ohll2006D will be seen as redshifted by 
the observer. 

However, while the integrated spectrum is certainly con- 
sistent with the presence of an outflowing bubble, it gives 
no indication of its physical size, or whether the bubble ge- 
ometry is appropriate. For example, the integrated spectra 
could be equally well explained by a large scale outflow, or 
a small bubble only just surrounding the starburst region. 
To investigate the energy of the outflow, and its ability to 
escape the host galaxy's gravitational potential, we must use 
spatial information in our spectra. 

We begin by considering the dominant contribution to 
the Lya line from photons travelling directly towards the 
observer. The wavelength of the blue cut-off in the line is 
determined by the outflow velocity seen by the observer. If 
the bubble is spherical and close to the galaxy we would ex- 
pect to see significant variation in the blue cut-off across the 
emission regio n since the li ne of sight velocity of the outflow 
depends on y/l — (b/R s ) 2 , where b is the impact parame- 
ter and R 3 is the shell radius. Clearly the consistency of 
the blue cut-off in Lya across the galaxy implies that the 
shell is much larger than the galaxy so that the illuminated 
portion of the absorbing screen has a planar geometry. We 
ill ustrate the likely view ing geometry in Fig. [7](cf., Fig. 18 
of IMas-Hesse et a! J 120031 ) for an outflow in the shell recom- 
binatio n phase (phase 4 in their classification scheme). In 
Haro 2. lLegrand et al.l l| 19971 ) are able to estimate the size of 
the shell from its low surface brightness Ha emission. They 
estimate that the projected diameter of the shell is 2.5 kpc, 
around five times the diameter of the emission region. A 
similar ratio would seem appropriate for the arc, although 
the shell could obviously be larger. 

Our data also require an additional contribution from 
photons scattered or created in the receding shell. For a large 
spherical shell, we would expect to see a widespread low 
surface brightness emission at this wavelength (e.g. Fig. [7]). 
However, no such distributed emission is seen in our data, or 
in Mass-Hesse et al.'s observations of local systems. Instead, 
by collapsing the datacube over a region dominated by the 
excess flux in the red wing, we find that the redshifted wing 
seems to have a similar angular distribution to that of the 
Lya emission coming directly from the galaxy. However, we 
can extend the model of Mass-Hesse et al. to explain the 
limited spatial extent of this emission as a result of the bi- 
conical outflow geometry. Redshifted photons from most of 
the receding shell must pass through the area of the galaxy 
that has not yet been ionised by the starburst. Even though 
they are redshifted by +u s hcii relative to the galaxy, they will 
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still be absorbed since the absorbing column is likely to be 
extremely high. Only in the ionised region surrounding the 
starburst will the Lya opacity be sufficiently low to allow 
photons to pass though the otherwise neutral disk (photons 
that pass through the ionised region are unlikely to be ab- 
sorbed by the foreground shell since they are redshifted by 
+2v s heii)- The column density of the residual gas disk is also 
likely to be much larger than that of the expanding shell. 
The result of this geometry is that the redshift wing will 
have similar spatial extent to the main Lya line. 

The proposed geometry explains the spatial distribution 
of the spectra without needing to invoke a special viewing 
angle. The only constraint is that the system is viewed down 
through the outflow cone. As the inclination between the 
outflow axis and the viewer becomes significant, we might 
expected to detect some rotation from the galaxy (a small 
shear is suggested by our SINFONI spectra of the [On] emis- 
sion), however, a shift of a several 10 kms -1 in the intrinsic 
profile of the Lya line is unlikely to be observable since the 
line shape is so strongly shaped by the blue-shifted absorp- 
tion (which does not share the velocity shear of the galaxy). 
If the system is viewed at a larger angle, Lya is likely to be 
completely absorbed since the expansion velocity perpen- 
dicular to the cone is small. In this case, it is unlikely that 
a redshift would have been obtained for this arc. It is also 
worth noting that the outflow spends much of its time in 
this post-blow-out phase. In earlier phases, the foreground 
shell will have lower column density or may even show emis- 
sion at its leading edge; in later phases, as the shell slows 
down and stalls, Lya does not escape at all. This model does 
not require us to be observing the galaxy at a particularly 
special phase of its evolution or at a special viewing angle. 

Whilst the proposed bi-polar outflow fits the observ- 
ables, there are other possible scenarios which would ex- 
plain the velocity offsets we measure: if the outflow were 
instead within the galaxy, and the Lya emission was non- 
Gaussian (e.g. a Lorentzian) then the red wing of emission 
seen at z =4.8860 could be explained without a receding 
shell. In this model, the blue-shifted material could arise due 
to an ISM outflow which is in the process of breaking out 
of the galaxy. However, the constancy of the Lya emission 
line shape and centroid is difficult to explain - any velocity 
motion within the star-forming regions should be reflected 
by the Lya, which is not seen in our data. 

4-1.1 Energetics of the outflow 

Although there are uncertainties in the geometry and col- 
umn density of the outflowing shell, it is nevertheless in- 
teresting to consider the energetics of the outflow, in or- 
der to see if the bi-polar outflow scenario is reasonable. For 
a normal stellar initial mass function, supernovae provid e 
~10 49 erg per solar mass of stars (eg.. iBenson et alj |2003). 
thus for a galaxy with a star- formation rate of 12±2M0 yr~\ 
~10 50 erg are available per year and the maximum extent of 
the Lya emission is ~ 12 kpc diameter. We estimate the 
distance of the shell from the galaxy by assuming that the 
linear size of the region covered by the swept up shell is at 
least q times the size of the maximum extent of Lya emis- 
sion region. For Haro 2, the shell is at least three times the 
spatial extent of the star- forming regions. If we adopt q < 3 
then we would expect to detect variations in the velocity of 



foreground absorber. This suggests a linear size for the end 
of the cone of £ 30 kpc. Adopting an opening angle for the 
outflow of ~ 60° , we estimate that the shell must be located 
Z 30 kpc from the starburst region. 

At a velocity of ~ 500 kms -1 it takes 60Myr to travel 
out to a distance of 30 kpc; within this time, supernovae will 
provide a total energy of ~6xl0 57 erg. Clearly there are a 
number of uncertainties in this estimate: not least that the 
shell may have decelerated from an initial higher velocity 
which would reduce the timescale. Nevertheless, we can es- 
timate the mass of the outflow via: M ou t flow = Nhi X Mh /x 
where Nhi is the column density, A is the area of the cone 
and Mh is the atomic mass of Hydrogen. Adopting the 
value of the observed column density we derived in i|3.2l 
of JVhi = 1.6 x 10 19 cm -2 and that the cone is uniform 
over an area of 700 kpc 2 , the total mass of the outflow is 
1.8 x 10 8 /a; Mq, where x is its Hi fraction. The kinetic en- 
ergy of the outflow is then E — \mv 2 — 5 x 10 56 fx erg. Thus, 
the outflow is energetically feasible if the neutral fraction is 
greater than 10%, or the column density of the shell is lower 
than we have assumed. (Note, however, that these two fac- 
tors well tend to play off against each other: a lower column 
density will tend to have a lower neutral fraction). 

It is interesting to note that the implied mass outflow 
rate is therefore ~ 3/x Mq yr~ . If we assume x=0.1, then 
the mass outflow rate is very efficient at ejecting baryons 
from the galaxy, with the mass loading of the wind being 
more than twice the star-formation rate of the galaxy. 

The estimates above are clearly uncertain. Nevertheless, 
it is interesting to see how far the expanding material might 
be expected to travel after the end of the starburst. We as- 
sume a halo circular velocity of 150 kms -1 and a wind speed 
of 500 kms -1 located 30 kpc from the center of the halo. If 
the wind flows freely out of the galaxy it is completely un- 
bound. However, the wind will slow down dramatically if 
it sweeps up material. To make an estimate of how far the 
wind will travel we assume that the wind expands into an 
NFW halo with a baryon fraction of 18% and that the wind 
has an opening angle of 60°. Even though the wind rapidly 
gains mass as it expands, it reaches a minimum of 160 kpc 
(almost 1 co-moving Mpc) before stalling. Thus, this star- 
burst driven wind has enough energy to pollute a volume of 
approximately ~3Mpc 3 . 



5 CONCLUSIONS 

The issues of metal ejection and feedback are some of the 
most outstanding problems in galaxy formation. One of 
the most important recent observational breakthroughs is 
that most high-redshift proto-galaxies appear to be sur- 
rounded by "superwinds" which are expelling material from 
the galaxy disc. However, at high redshift (where star forma- 
tion efficiency was at its peak and therefore this phenomenon 
is likely at its peak activity), observational evidence is usu- 
ally based on measuring the velocity offsets between the neb- 
ular emission (such as [On] or Ha) and rest-frame UV lines 
such as Lya or UV ISM lines. These observations are usually 
based on longslit observations which lack spatial informa- 
tion, which is crucial if we are to understand the dynamics 
and fate of the outflowing gas. 

In this paper, we have spatially resolved and mapped 
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Figure 7. This figure illustrates the details of the proposed configuration to explain the spatially resolved emission spectrum. The Lya 
emission comes predominantly from photons generated in the star forming region 'A'. The Lya emission region is physically larger than 
the due to recombination in dense gas clumps that have not yet been swept up and ionised by the outflow, and the resonant Lya may not 
escape readily from the central regions of the starburst. The outflow cavity 'B' contains predominantly ionised material flowing out of 
the galaxy. It is hot and diffuse and therefore not seen directly. However, the shell 'C (travelling towards the observer) that has formed 
as the wind has swept up the interior of the cavity can be seen in absorption at — u a hcll against the central Lya emission region. The 
geometry of the outflow, and its extent compared the source mean that the projected velocity is almost constant across the illuminated 
portion of the shell. Photons that scatter of the receding shell ('D'), or are produced by ionisation of its inner edge by UV photons 
from the starburst, are seen by the observer as a weak emission line (broadened by radiative transfer effects) at +u a hell- The bi-conical 
nature of the outflow is key to explaining why the redshifted emission is only seen over the portion of the shell directly behind region 
'A'. Redshifted photons from other regions of the receding shell must pass through the area of the galaxy that has not yet been ionised. 
Even though they are redshifted relative to the galaxy, they will be absorbed since the absorbing column is likely to be extremely high. 
Note that photons that pass through the central ionised region of the galaxy are unlikely to be absorbed by the foreground shell since 
they are redshifted by +2u a h e ii- 



the dynamical properties of a highly magnified gravitation- 
ally lensed galaxy at z=4.88. The main features are de- 
scribed below: 

• The lensing amplification is a factor 16±2 
(Am=3.0±0.2 magnitudes). In the (reconstructed) source- 
frame the HST imaging suggests the galaxy is only ~ 
2 x 0.8 kpc FWHM. Moreover, the lensing corrected /-band 
magnitude is 7=25.2. For comparison, an L* Lyman-break 
galaxy at z ~ has i=25.4, suggesting this galaxy is typical 
of galaxies at these early times. Due to the lensing ampli- 
fication we are able to study a typical z ~5 galaxy at high 
signal-to-noise and spatially resolve the galaxy on <;200pc 
scales. 

• The [Oll]AA3726.2,3728.9 emission line maps the (sys- 
temic) nebular emission and shows a peak-to-peak velocity 
gradient of ^60kms _1 across 3 kpc in projection. The inte- 
grated [On] emission line flux suggests an integrated star- 
formation rate of 12±2Mq yr _1 . The velocity shear, spatial 
extent and velocity dispersion of this nebular emission sug- 
gest a dynamical mass of ~ 1 x 10 10 Mq within a radius of 
2 kpc. 

• The Lya emission is much more extended than ei- 
ther the UV continuum or [On] emission, extending over 
11.9x2.4 kpc (FWHM). The Lya emission has a asymmet- 



ric profile causing it to appear redshifted from the systemic 
velocity by +200±40kms~ 1 . Across the whole galaxy im- 
age both the shape and centroid of the Lya emission are 
remarkably constant with the centroid of Lya varying by 
<10kms _1 . 

• Both the collapsed and spatially resolved Lya emis- 
sion are best fit with a 3 component model: (i) an underlying 
(Gaussian) emission profile which has a redshift and width 
in excellent agreement with the nebular emission; (ii) an 
absorber with a Voigt profile at — 500 km s -1 and a column 
density of 1.6x 10 19 cm 2 , and (iii) a broad redshifted emission 
line component at +400 km s -1 . 

• In the collapsed rest-frame UV spectrum, the weak 
UV ISM lines of Siiv are blue-shifted from the systemic ve- 
locity by — 400TT00 kms -1 (the velocity of these lines are 
in agreement with the velocity of the absorber seen against 
Lya). 

The velocity offsets between the nebular emission and 
spatially extended Lya are comparable to those seen in most 
high redshift LB Gs at similar redshifts (|Shaplev et al.ll2003l : 
lErb et al.ll2003l) and are usually attributed to galactic-scale 
outflows. Hence, this leads us to interpret the results in the 
context of an evolving starburst outflow and we illustrate 
the likely viewing geometry in Fig. [7] In this model, the [On] 
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emission is assumed to trace the underlying (star-forming) 
region embedded in an bi-conical outflowing bubble in which 
the starburst generates an over-pressurised region within the 
ISM which becomes Rayleigh- Taylor unstable in which the 
initial shell breaks up, leaving behind denser clumps which 
may be responsible for the extended Lya region. Any outflow 
will break out of the galaxy along the axis of least resistance 
and hence for gas in a planar or bar-like galaxy this generates 
an approximately bi-conical outflow. 

The proposed geometry naturally explains the spatial 
distribution of the spectra and the velocity offsets and struc- 
tures we observe. We consider the energetics of the outflow 
and show that the proposed model is self-consistent and en- 
ergetically feasible. The energy of the wind is such that it 
plausibly reaches ~160kpc (almost 1 co-moving Mpc) be- 
fore stalling. At this point, it will have polluted a comoving 
volume of nearly 3 Mpc 3 with metal enriched material. 

Whilst these observations are based on a single galaxy, 
our results clearly show the power of combining gravitational 
lensing with optical and near-infrared observations to probe 
the star-formation activity, masses and feedback processes 
in typical high redshift galaxies in great detail. The next 
step is to generate a statistically useful sample to gauge 
the prevalence outflows from these young galaxies. In turn 
these measurements will help us understand the wide-spread 
enrichment of the early universe, and may explain why only 
10% of baryons cool to form stars. 
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APPENDIX A: SERENDIPITOUS 
BACKGROUND GALAXIES 

Since we are able to simultaneously survey all of the critical 
lines from 2=0 to 2=7 with the VIMOS IFU and part of the 
2=1.5 to 2=14 critical lines with the SINFONI IFU we ex- 
ploit the SINFONI and VIMOS data to search for serendip- 
itous sources behind the lensing cluster. We identify seven 
candidates between z =1.5 and z =5.5 from both optical 
and near-infrared spectroscopy, at least three of which have 
either optical or near-infrared broad-band counterparts. One 
of these sources is a radial counter-image of the 2=4.88 arc 
which is then used to provide strong constraints on the lens 
modelling (j) l3.1|) . The remaining sources in the VIMOS dat- 
acube are identified as [On] emission at 2=0.99 and Lya at 
2=3.66. The SINFONI IFU emitters are tentatively identi- 
fied as either Ha at 2=1.5 to 2=2.7 or [On] at z ~5.3 (Ta- 
ble [XT}. These observations show the power which integral 



field spectroscopy will have in the future in finding high- 
redshift galaxies b ehind cluster cores (e.g. lEllis et al.ll200ll ; 
iKneib et al.ll2004T ). The observations presented here also al- 
low an insight of the deeper observations that will be made 
possible with even larger area integral field spectrograph's 
on th e VLT in the near future (e.g. MUSE; | Henault et al.l 
120041) or the mosaic mode of KMOS l|Sharples et alj|2004h '). 
These observations will open the doorway to probing (and 
spatially resolving) the properties of the primeval galaxies 
responsible for reionization around 2=8-12. 



Al Identifying Serendipitous Galaxies in the 
datacube 

Using the IFU coverage of the cluster cores, we attempt 
to identify serendipitous sources in the IFU datacubes. For 
the wavelength range of VIMOS (5000A-l/xm), we are able 
to search for Lya between 2=3.11 and 2=6.8 and [On] be- 
tween 2=0.34 and 2=1.57. The wavelength coverage of the 
SINFONI IFU observations (A=1.451-2.463^m) also allow 
us to search for [On] emission between redshifts z =2.90 to 
2 =5.61, Ha between z =1.21 and z =2.75 and even Lya 
between z =10.9 and z =19.3). 

To identify serendipitous line emission we scan each 
pixel of the datacube. At each wavelength in each spectrum 
we calculate the noise over ±150A and demand a signal-to- 
noise of four to detect an emission line in at least three of 
the eight adjacent pixels (we mask out any regions of strong 
sky emission from our analysis completely) . If line emission 
is detected, we then fit the emission line with a single Gaus- 
sian profile and require that the emission line must have a 
line width greater than that of the sky at the given wave- 
length. Using this technique, we identify three candidates 
in the VIMOS IFU datacube (excluding [O n], [Om] or H/3 
emission from the cluster members) and four candidates in 
the near-infrared (SINFONI) datacube. We show the one 
and two-dimensional spectra (as well as their location in the 
datacube) in Fig. lAlfclA2l and briefly review their properties 
and likely redshifts here. 
Optical 

VIMOS Source 1 (Rl): The first serendipitous source is 
the radial counter-image of the 2=4.879 arc (labelled Rl 
in Fig. [JJ . The position of this counter-image is used as an 
extra constraint in the lens modelling in § 13.11 

VIMOS Source 2 (Al): In the HST imaging of the clus- 
ter, the third counter-image of the 2=4.789 arc (at approx- 
imately +6", +14") shows a bright knot superimposed onto 
the arc which does not appear in any of th e other counter- 
imag es of the arc and therefore (as noted bv lGladders et all 
2002) is most likely a foreground galaxy super-imposed 
onto the arc. This galaxy has a U-band magnitude of 
V606=25.6±0.2 and we identify an emission line at 5665A 
with an emission line flux of 1.4±0.5x 10 _16 erg s _1 cm~ 2 and 
an observed equivalent width of 700±150A. If the line is 
identified as [On], then a redshift of 0.51 yields an EW of 
~470A, which is high compared to local star-forming galax- 
ies. We therefore suggest that the emission is Lya at a red- 
shift of 2=3.66 with a rest-frame EW of 150±50A. The emis- 
sion line width (FWHM=12.5A) suggests a line width of 
220±60 kms" 1 . 

VIMOS Source 3 (A 2): The third background galaxy de- 
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Figure Al. Spectra of the three background, serendipitous 
sources in the VIMOS IFU field of view. In all three pan- 
els we show the collapsed, one-dimensional spectrum, the two- 
dimensional spectrum (generated by unwrapping the datacube 
around the object) as well as the location of the object in the 
datacube (right). Top: The super-imposed foreground galaxy (la- 
belled Al in Fig. 1) on the third image of the 2=4.88 arc has 
strong line emission at 5600A which we identify as Lya yielding 
a redshift 2=3.66; Middle: The 2=4.789 arc is a radial counter- 
image (labelled Rl in Fig. 1) of the arc and is used as a primary 
constraint on the lens model in § 13.11 Bottom: The blue-disk 
galaxy (labelled A2 in Fig. 1) is lensed source at 2=0.99 from 
identification of [On] at 7400A. 



tected is the blue disk-like galaxy labelled A2 in Fig. [T] We 
detect strong emission at 7409A but no other significant 
(>5cr) emission lines in the one-dimensional spectrum. The 
emission line also has a velocity shear of ~5A across the 
~1.5". If the emission is identified as [On]A3727 a redshift 
of 0.9855 is derived, placing the galaxy behind the cluster, 
and probably lensed. However, it is also possible that the 
emission is Ha at z =0.129, although we place limits on the 
emission line ratios of Ha/H/3>12 and Ha/[Oin]A5007>6. 
Identification as [Om]A5007 is ruled out due to the non- 
detection of [Oiii]A4959. Alternatively, if the emission is 
H/3 at z =0.524, we constrain the emission line flux ra- 
tios as H/3/[Om]A5007>8 and H/3/[0n]A3727>17. We there- 
fore suggest the most likely identification is [Oll]A3727 at 
z =0.9855. 
Near-Infrared 

SINFONI Source 1: The proximity of this source to the 
z=4.88 critical line makes it more likely to be [On] for 
z=5.25. We can not rule out Ha at 2=2.55, although the 
lack of [Oiii] or H/3 at 1.8^im (although near to the atmo- 
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Figure A2. Spectra of the four background, serendipitous 
sources in the SINFONI IFU field of view. In the left hand panels 
we show the one-dimensional spectra with the two-dimensional 
spectra below. In the right hand panel we show the two- 
dimensional (narrow-band) images which show the positions in 
the datacubes for these sources. The positions of these sources 
are also shown in Fig. 1 and details are also given in Tabic lATI 
The images in the right hand panel are generated by collapsing 
the datacubes over the wavelength range shown by the two dashed 
lines in the one-dimensional spectra. 



spheric Jf-band absorption) or [Nll]A6583 emission makes it 
more likely [On]. 

SINFONI Source 2: This source lies close to the high ampli- 
fication lines for z=2-3 and therefore we suggest Ha is the 
most likely line identification. There is no strong [Oiii] in 
the if-band. 

SINFONI Source 3: The proximity to the z~1.5 critical lines 
makes it more likely to be Ha than high redshift [On] 

SINFONI Source 4-' This is the weakest emission line in the 
sample, although it still fulfils the selection criterion. The 
most likely line identification is [On] at z=5.51 due to its 
proximity to the z ~5 critical curve. This emission line could 
also be identified at Ha for z=2.70, however, we note that 
there are no signs of strong [Oiii] with a flux ratio limit of 
[Om]/Ha emission line ratio £ 0.2 



APPENDIX B: SPECTROSCOPIC 
PROPERTIES OF THE CLUSTER GALAXIES 

Table [Bl gives the spectroscopic properties (RA, Dec, [On] 
emission line flux, line width (FWHM) and Equivalent 
Width (EW)) of the cluster galaxies. 
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Table Al. 

Properties of the Serendipitous Sources 



ID 


RA 


Dec 


^obs 


Line Flux 


Likely redshift 


observed 








(jim) 


(xlO-^ergs" 1 ) 




(A) 


VIMOS 














Al 


02:24:33.862 


-00:02:17.61 


0.5665 


1.4±0.5 


3.660 (Lya) 


470±60 


A2 


02:24:34.602 


-00:02:34.80 


0.7400 


6.0±1.0 


0.986 ([On]) 


102±20 


Rl 


02:24:34.466 


-00:02:34.77 


0.7148 


4.5±0.8 


4.879 (Lya) 


720±30 


SINFONI 














1 


02:24:33.455 


-00:02:21.73 


2.3330 


0.1±0.05 


2.554 (Ha) / 5.259 ([On]) 


>78 


2 


02:24:33.591 


-00:02:24.22 


2.2455 


0.06±0.03 


2.421 (Ha) / 5.024 ([On]) 


>47 


3 


02:24:33.708 


-00:02:25.71 


1.6740 


0.04±0.02 


1.550 (Ha) / 3.159 ([Oil]) 


>31 


4 


02:24:33.315 


-00:02:23.28 


2.4302 


0.03±0.02 


2.702 (Ha) / 5.520 ([Oil]) 


>23 



Table Al. Properties of the serendipitous sources in the VIMOS and SINFONI IFU fields of view 



Table Bl. 

Spectroscopic Properties of the Cluster Galaxies 



ID 


RA 


Dec 


2 


Features 


[On] Flux 


[On] FWHM 


EW 












(xl0" 17 crg/cm 2 /s)) 


(A) 


(A) 


CGI 


02:24:34.097 


-00:02:30.90 


0.7784 


H&K 








CG2 


02:24:34.246 


-00:02:32.54 


0.7725 


H&K 








CG3 


02:24:34.184 


-00:02:31.29 


0.7603 


[On] 


8.5±1.0 


5.7 


7.0 


CG4 


02:24:34.717 


-00:02:35.22 


0.7845 


[On] 


12.0±1.0 


7.0 


8.2 


CG5 


02:24:34.751 


-00:02:28.75 


0.7864 


H&K, [On] 


10.0±2.0 


4.1 


22.0 


CG6 


02:24:34.472 


-00:02:27.71 


0.7848 


H&K, [On] 


20.0±1.0 


11.1 


17.0 



Table Bl. Note: Properties of the brightest cluster galaxies, as measured from our VIMOS IFU observations. Equivalent widths (EW) 
and emission line FWHM are given in the rest-frame of the galaxy. 



